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Abstract 



We explore what can be said on the effective temperature and sound speed of a statistical 



> 

(N 

m 

00 

O ' ensemble of fluid phonons present at the onset of a conventional inflationary phase. The 

phonons are the actual normal modes of the gravitating and irrotational fluid that dominates 
the protoinflationary dynamics. The bounds on the tensor to scalar ratio result in a class 
of novel constraints involving the slow roll parameter, the sound speed of the phonons and 



H ■ the temperature of the plasma prior to the onset of inflation. If the current size of the 

" " Hubble radius coincides with the inflationary event horizon redshifted down to the present 

epoch, the sound speed of the phonons can be assessed from independent measurements of 

the tensor to scalar ratio and of the tensor spectral index. 
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There are some who argue that the primeval plasma went through a phase of quasi-de 
Sitter expansion whose likelihood is difficult to gauge given our ignorance of the preinflation- 
ary initial conditions. Since inflation cannot be anyway eternal in its past, the inflationary 
stage of expansion is customarily complemented by a preinfiationary epoch of decelerated 
evolution wherqj a > and a < 0. A decelerated phase driven by an irrotational fluid 
possesses its own normal modes that depend on the sound speed of the medium and on the 
specific coupling to the background geometry. The quantum theory of the fluctuations of 
gravitating, irrotational and relativistic fluids [1] has been developed even prior to the actual 
formulation of the conventional inflationary paradigm and in the context of the pioneering 
analyses of the relativistic theory of large-scale inhomogeneities j2]. The idea pursued here- 
under is, in short, to scrutinize the implications of an ensemble of protoinflationary phonons 
for the large-scale gravitational perturbations and to compare the obtained results with the 
standard lore provided by the concordance paradigm. 

The preinfiationary physics is heuristically modeled by often assuming that the inflaton 
is the one and only source of the geometry while the initial conditions of the mode func- 
tions are taken to differ from the vacuum |3] (see also ^). We ought ffist to relax the 
aforementioned assumption by admitting that initially the total energy-momentum tensor is 
dominated by an irrotational fluid so that the approximate forms of the scale factors across 
the protoinflationary transition can be presented in the form of a continuous interpolatioio: 
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,,{t) = a,{H,tr, t<H;\ 0<a<l, (1) 

ainf(t) = a*[(a/7)^*i + (7-«)/7]^ t > H;\ 7 > 1. (2) 

The time scale H~^ defines the transition between deceleration and acceleration but three 
subsequent regimes can be physically distinguished: ffist the preinfiationary epoch (where the 
fluid dominate); second the protoinflationary regime (where a changes sign); third the fully 
developed inflationary stage. Defining e = —H/H"^ as the rate of variation of the Hubble 
parameter, we shall have that e c^ 3{w + l)/2 > 1 during the ffist regime. Conversely 
e ~ 0{1) during the protoinflationary phase and e ^ 1 as soon as the slow-roll evolution 
starts. Explicit examples of this dynamics can be found in terms of either analytic or 
numerical solutions of the Einstein equations with multiple sources. Sticking to a simplistic 
model, the interpolating scale factor a{t) = a*(sinh/3if*)^/^ (with /3 = 3{w + l)/2) is a 



^In the present analysis the cosmic or the conformal time parametrizations shall be used indifferently. The 
derivative with respect to the cosmic time t will be denoted by the overdot, while the derivation with respect 
to the conformal time coordinate t will be denoted by a prime. The relation between the parametrizations is 
given by dt = a{T) dr where a{T) is the scale factor of a Friedmann- Robertson- Walker background geometry 
in the conformal time parametrization. Consequently, the Hubble rate H — (d/a) coincides with T-l/a 
where T-L = {a' /a). Spatially flat background geometry shall be mainly discussed to make contact with the 
concordance paradigm but the obtained results can be easily extended to the spatially curved geometries. 

•^If the protoinflationary phase is dominated by an irrotational fluid with constant barotropic index w, 
a = 2/[3(w + f )]. Note that the interpolation dl])-© is continuous in t, = H^^ , i.e. aproto(i*) ~ ainf(i*) 
and apioto(4) = ainf(4). 



solution of the Friedmann equations supplemented by a perfect barotropic fluid and and by 
a minimally coupled scalar field with inflating potential; in this case e(t) = /3/cosh^ {f3H^t). 
During the protoinflationary phase the fluctuations of an irrotational and perfect fluid 
are described in terms of an action whose explicit form has been derived, for the first time, 
by Lukash [T] 
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where 7*-' = 5^^ in the spatially flat case and clij) = p[/p[ denotes the total sound speed of 
the system measured in natural units. The variable gph of Ref. [I] can be expressed in terms 
of 7^, i.e. the curvature perturbations on comoving orthogonal hypersurfaces satisfying, 
because of the action ([3]), the following equation: 



^// ^ 2^i^7^' - cj v^7^ = 0, n = ~qph/z, 
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Besides the fluid phonons the other modes of vibrations of the protoinflationary fluid 
are associated with spin-two inhomogeneities [2] whose action follows from the second-order 
fluctuation of the Einstein-Hilbert action in the tensor amplitude hij (with dik^^ = h] = 0), 
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where £p denotes the Planck length, Mp is the Planck mass and Mp is the reduced Planck 
mass. The variables hij, TZ and gph are gauge-invariant and can be evaluated in the coordinate 
system which is more suitable for the problem at hand. This statement holds, of course, for 
the scalar modes since the tensors are automatically gauge- invariant. Using, for instance, 
the off diagonal gauge [5] (where (5s5'oo = 2a^0 and 45'oi = —a^P) the fluctuations of the full 
metric are: 
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In Eqs. (jSD^dZD the tensor and scalar perturbations are denoted, respectively, by 6t and 6s; 
6sgij{x,T) vanishes in the off diagonal gauge [5j. Equations (jHD follow from the perturbed 
Einstein equations implying, in the absence of anisotropic stresses, (3' + 2'Hf3 + = 0; the 
Bardeen potential |6], in the off-diagonal gauge, can be expressed as \E'(x, r) = — 'H(r)/3(af, r). 
After promoting the normal mode qph{x, t) and hij{x, r) to the status of field operators 
obeying (equal time) canonical commutation relations with their conjugate momenta, the 



mode expansions of the corresponding operators become: 



hij{x,T) 

'J 1 _ 



(27r)3/2 



fkir)a^,e-^'-^ + f:iT)dU^'-^ 






^(k) 



K\9k,\{r)e 



-ik 



''' + b\x9UVe 



ik-x 



(8) 
(9) 



where [%_;^, 6J-,a'] = ^^'^K^ - p}5xx' and [d^, at] = 5'^^'^{k 



p); eij (k) accounts for the two 



propagating polarizations of the tensor modes and satisfies e^ e^ 



26 \y. The mode 

functions fk and gk,\ obey two similar equations leading, however, to different dispersion 
relations: 
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During the protoinfiationary regime e decreases from values e(t) > 1 to values e(t) ~ 0{1). 
As soon as e(t) ^ 1, the evolution of curvature perturbations is governed virtually by Eq. 
(TTTj) except for the dispersion relations (implying Cg — )■ 1) and except for a different form of 
the pump field i.ejj 



n" + 2^7^' - v^n = o, 
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where ip denotes the infiaton field and q^ is just the analog of gph of Eq. ([3]). As discussed 
after Eqs. ([I])-(l2l) a smooth protoinfiationary transition implies the smoothness of the scale 
factors, of the curvature pump fields (i.e.2;ph and z^) and, ultimately, of the extrinsic cur- 
vature. The inhomogeneities can be matched by requiring the continuity of the induced 
three metric and of the extrinsic curvature on the space-like hypersurface of the transition 
[8]. Such a requirement demands, at the level of Eqs. ([6])-([7]), the continuity of 71{x,t), 
/3{x, t) and of hij{x, r). The results obtained in this way can be directly checked by deriving 
the evolution equations during a phase where the fluid phonons and the infiaton quanta are 
simultaneously present but the fluid component is already subleading (see e.g. last paper of 
Ref. [3]). 

Since the phonons and the gravitons evolve, initially, in a fluid background (which may 
even be dominated by radiation) the initial state for each species must be reasonably de- 
scribed by a mixed density matrix with Bose-Einstein distribution i.e. 
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where n^ is the average multiplicity of each Fourier mode and \{n}) = |n^^) ® jn^^) ® l^fc )••• 
with the ellipses standing for all the occupied modes of the fleld. The average multiplicities 



"^Equation ()lip can be derived within the same logic employed in the first paper of Ref. W but in the 
case of scalar field matter [7] . 



of gravitons and phonons are given by: 
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having denoted with k^ = T the effective temperature of the phonons and of the gravitonqj. 
Note that u^^ = kcg while u^'^ = k. Muhiplying the occupation numbers of Eq. ( TT3l) by 
the corresponding quanta of energy and integrating the result over the appropriate phase 
space, the energy density of phonons and gravitons becomes ges'^^T^/'^O where g^s denotes 
the effective number of relativistic degrees of freedom. If only gravitons and phonons are 
considered ges = 3 but g^s maybe larger than 3 because of other species in local thermal 
equilibrium. If Tmax = [45/(47r^(yfefF)]^/'^-y/if*Mp the energy density associated with the initial 
state is 0{3H^Mp). For the forthcoming apphcations it will be practical to define T^, = 
QT^ax where Q < I measures the protoinflationary temperature in units of Tmax- 

The initial mixed quantum state of the system is the direct product of the quantum state 
of the phonons and of the gravitons, i.e. |$) = iV^phonons) ® |0gravitons) and the corresponding 
power spectra are obtained from the following two expectation values: 

($|7^(f,r)7^(f +f,r)|<|.) = f ^Vn{k,rf-^, (14) 

J k kr 

^ ^ [' H jc Sin nf 

($|/i,,(f,r)/i^^(f + f,r)|$) = y — PT(A:,r)-^, (15) 

wherqj the two power spectra are: 

VAk,r) = —,\fk{rn2nt + l), Vr{k,r) = -^^ k' \g,{r)\\2nf + 1) . (16) 

These spectra must be evaluated in the limits k/{a^:H^) ^ 1 (corresponding to wavelengths 
smaller then the Hubble radius at the protoinflationary transition) and k/{aH) ^ 1 (since 
the phenomenologically relevant wavelengths are still larger than the Hubble radius at the 
epoch of the matter radiation equality). In these concurrent limits Eq. flT6l) implies 

8 /V\ / k \"s~i / kc \ 

Vnik,r) = ^[j)[^) -F^(fc,r.,c.,e,r,)coth(^^j, (17) 

128 /V\/fc\"T f k \ 

^-<*-^' = -(m?)(^) ^T(*.r..c„.)coth(— ). (18) 

where rj = (p/{H(p) = (e — r/) and V is the infiaton potential; the two spectral indices are: 

ris = 1 - 6e + 2r7, nx = -2e. (19) 



■'We shall posit that phonons and gravitons have the same temperature but this assumption can be 
dropped if only approximate (kinetic) equilibrium holds between the different species. 

^We used, in Eq. p3|) that the mode function of each tensor polarization is the same since the two 
polarizations do not interact between them. 



For k/{a^H^) ^ 1 the two functions J-'-ji and J-t go to 1 up to corrections 0{k~'^T~'^) which 
will be irrelevant for the subsequent discussion. From Eqs. ( 1T7|) and (TTS!) the tensor to scalar 
ratio is given by 

where the power spectra have been referred to the standard pivot scale fcp = 0.002 Mpc~"^ 
corresponding to an effective multipole fgff — 30. 

In the concordance paradigm radiation takes places almost suddenly after inflation and 
long phases (different from radiation) are excluded down to the curvature scale of matter- 
radiation equality. Under this approximation the ratio (/c/Zct) is given by 



where iV^ax denotes the maximal number of efolds presently accessible through large-scale 
observations and Ar is the amplitude of the power spectrum of curvature perturbations 
determined at the pivot scale kp. In practice A^max is determined by fitting inside the Hubble 
radius the approximate size of the inflationary event horizon redshifted down to the present 
epoch: 

iVmax = 62.2 + 0.5 In (10^ - In (ho/O.T) + 0.25 In [10^ /i^nRo/4.15], (22) 

where ^ = neAn and Qjiq is the present critical fraction of radiation energy density. If 
the postinfiationary phase contains long epochs expanding with a rate which is slower than 
radiation A^max can even be 0{15) efolds larger. 

The WMAP9 yr data (first and second papers of Ref. J9]) alone imply r'Yikp) < 0.38. 
When the WMAP9 data are combined with other data setCl we have that the cosmological 
parameters are determined to be 

(l^bo, ^cc^deo, ho, ris, e^e) = (0.04596, 0.2360, 0.7180, 0.6969, 0.9647, 0.081), (23) 

while A-ji = 2.431 X 10^^. In the case of the parameters of Eq. (12 3 p the bounds on r^ and on 
the tensor spectral index imply, respectively, VT^kp) < 0.13 (95 % CL) and tt-t > —0.016. The 
Planck satellite gives r^ikp) < 0.11 (see the last paper of Ref. [10]) which is quantitatively 
comparable with the WMAP9 limit. 

If A^ > A/'max Eqs. ([20]) and (^Hj-I^ imply that r^ikp) -^ 16e. Similarly, for c, ^ 1 (i.e. 
when the sound speed coincides with the speed of light) rT:{kp) = 16e in spite of the total 
number of efolds. This is the so-called consistency regime stipulating that, at the pivot scale, 
tt — 16e ~ — 8?2t- Whenever the speed of sound does not equal the speed of light the consis- 
tency relation is violated. According to Eqs. (12T]) -( 123]) . if A^ ~ C'(A'max) then k/k^ ~ C(l) as 



^In particular the data on the Hubble rate, the ones on the baryon acoustic oscillations supplemented 
by the data of the Atacama Cosmology Telescope, by the data of the south pole telescope and by the three 
year sample of the supernova legacy survey (see the papers reported in Ref. 10: ) 
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Figure 1: The tensor to scalar ratio for the fiducial set of parameters of Eq. (123!) . In the left 
plot Q = 0.01 and k is allowed to vary; in the right plot k = k^. On the vertical and the 
horizontal axes the common logarithms of the corresponding quantities are illustrated. 



long as 60 < A^ < 64, 0.01 < Q < 1 and 0.26^^^^^ < (fep/fcx) < 0.48^^^^^ Thus r^ikp, e, c^, A;t) 
can be expanded for kp/{2k-[) < 1 so that, approximately, rx — 16eCs(Tgi./Tph) where the 
case of different (kinetic) temperatures has been included for completeness (in the thermal 
equilibrium case Tph — T^,.). This means that the consistency relation is violated and the 
sound speed of the primordial phonons can be potentially measured from independent esti- 
mates of the tensor to scalar ratio and of the tensor spectral index since it is still true that 
riT = — 2e. In Fig. [T]the value of tt is illustrated without approximations and for different 
values of the inflationary efolds. In the left plot Q is held flxed while in the right plot k co- 
incides with fcp. In both plots the sound speed has been chosen to equal the radiation value, 
i.e. Cg = 1/3. With the squares we illustrate the WMAP9 limit on rx. The allowed region 
is always below the squares. The WMAP9 limit is purely illustrative since it is derived, as 
usual, by assuming the consistency relation. 

All in all the results of this investigation imply that a protoinflationary phase dominated 
by fluid phonons invalidates the consistency paradigm by impacting on the tensor to scalar 
ratio whose spectral dependence is determined by the sound speed of the phonons and by the 
potentially different temperatures of the various species in kinetic equilibrium. Reversing the 
argument, the present flndings show more modestly that the consistency relations, even if 
heuristically assumed in nearly all experimental analyses, are not always tenable and bear the 
mark of an essential completion of the early inflationary dynamics, i.e. the protoinflationary 
transition. It is therefore tempting to speculate that independent measurements of tt and 
riT from B-niode polarization may offer a potentially novel diagnostic of the role played by 
primordial phonons in setting the initial conditions of large-scale gravitational perturbations. 
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